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Knowledge of the binary population in stellar groupings provides important information about the outcome of the star 
forming process in different environments. Binarity is also a key ingredient in stellar population studies and is a prerequisite 
to calibrate the binary evolution channels. In these proceedings we present an overview of several commonly used methods 
to pair individual stars into binary systems, which we refer to as the pairing function. Many pairing functions are frequently 
used by observers and computational astronomers, either for the mathematical convenience, or because they roughly 
describe the expected outcome of the star forming process. We discuss the consequences of each pairing function for the 
interpretation of observations and numerical simulations. The binary fraction and mass ratio distribution generally depend 
strongly on the selection of the range in primary spectral type in a sample. These quantities, when derived from a binary 
survey with a mass-limited sample of target stars, are thus not representative for the population as a whole. 
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1 Introduction 

Most stars form as part of a binary or multiple system (e.g., 
Duquennoy & Mayor 1991; Mason et al. 1998; Goodwin & 
Kroupa 2005; Kouwenhoven et al. 2005, 2007b; Goodwin 
et al. 2007), and about 15% of these are part of a multiple 
{N > 3) system (Tokovinin & Smekhov 2002; CoiTeia et 
al. 2006; Tokovinin et al. 2006; Hu et al. 2008). In order to 
understand the star formation process, it is of crucial impor- 
tance to characterize the properties of the binary population. 
In these proceeding we discuss several methods of pairing 
stars into binary systems, and refer to these as pairing func- 
tions. The pairing functions described in these proceedings 
are used in literature for (i) a description of the observed 
binary population, (ii) initial conditions for iV-body sim- 
ulations, and (iii) a description of the outcome of hydro- 
dynamical simulations or theoretical predictions of the star 
formation process. The pairing function determines the dis- 
tribution of stellar masses. Variations of the primordial pair- 
ing function with environment thus imply variations of the 
initial mass function (IMF) and the star formation process. 

2 Method and terminology 

Consider a stellar population consisting of single stars and 
binary systems. The pairing algorithm used to generate the 
binaries in such a stellar population is called the pairing 
function (Kouwenhoven et al. 2008). This pairing function 
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could, for example, be random pairing of the binary compo- 
nents from the initial mass function (IMF). 

For a binary system consisting of two stars with masses 
Ml and A/2 < Mi, we refer to Mi as the primary star and 
M2 as the companion star. The mass ratio of the system is 
q = M2/M1 with < q < I, and the total mass of the 
system is Mt — Mi + M2. For simplicity, we neglect the 
presence of triple and higher-order systems. 

Depending on the pairing function, the mass of the pri- 
mary star or the total mass of each binary system is drawn 
from a generating mass distribution fniM), such as the 
Kroupa (2001) mass distribution. A certain fraction B of 
the stars is assigned a companion star. Depending on the 
pairing function, additional manipulations may performed, 
such as swapping of primary and companion stars such that 
Ml > M2, or the rejection of low-mass companions (see 
§|3j- Note that, due to these manipulations, the stellar mass 
distribution (of all stars: primary, companion, and single 
stars) changes. The resulting mass distribution ("IMF") is 
therefore generally not equal to fniM). The only excep- 
tion is random pairing (RP), where both components of each 
binary system are independently drawn from f^i{M). 

3 Several pairing functions 

There are numerous ways to pair individual stars into bina- 
ries. Several straightforward pairing functions are discussed 
below (we refer to Kouwenhoven et al. (2008) for details). 
- Random pairing (RP). Both components of each binary 
are drawn from a generating mass distribution fuiM) and 
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paired into a binary system, and swapped, if necessary, so 
that the primary is the most massive of the two (e.g., Piskunov 
& Mal'Kov 1991; Mal'Kov & Zinnecker 2001). RP is the 
simplest possible pairing function, as it is fully characterised 
by /m (M) and B. RP is the only pairing function for which 
the resulting stellar mass distribution fM.ai\{^I) of all stars 
is equal to fmiM). Random pairing in a stellar population 
implies a chaotic formation process, where stars are paired 
into binaries irrespective of their mass. 

- Primary constrained random pairing (PCRP). The pri- 
mary mass Ml is drawn from fM{M). The companion star 
is also drawn from fM{M), with the additional constraint 
that it is less massive than the primary: M2 < Mi. Note the 
subtle difference between RP and PCRP. 

- Primary-constrained pairing (PCP). The primary mass 
Ml is drawn from the mass distribution fuiM), and the 
mass ratio q is drawn from a mass ratio distribution fq{q). 
Finally, the companion mass M2 — qMi is calculated. A 
possible origin of this pairing function comes from disk 
fragmentation. The mass of a circumstellar disk is related 
to the mass of the primary star 

- Split-core pairing (SCP). The total mass Mr = M1+M2 
of the binary system is drawn from a generating mass dis- 
tribution, and the mass ratio of the binary system is drawn 
from a mass ratio distribution. Subsequently, Mt is split up 
into two stars with a mass Mi = Mt(1 + q)^^ for the 
primary and a mass M2 = Mt{1 + q~^)~^ for the compan- 
ion star The resulting binary system than has a total mass 
Mt = Ml + M2 and a mass ratio q — M2/M1. A pair- 
ing function similar to SCP is expected from the theory of 
fragmenting gas clumps. 

For pairing functions PCP and SCP an additional com- 
plication may occur. If a low mass ratio q is drawn from 
fq{q), it is possible that the companion mass is smaller than 
a certain minimum mass M^in- Such a minimum mass may 
originate from star formation theories. A minimum compan- 
ion mass also originates from our definition of a binary sys- 
tem: a star with a companion less massive than the deuterium- 
burning limit (i.e., a planet) is considered a single star, im- 
plying Mniin ~ 0.02 M0. Adopting this lower limit would 
also imply that stellar and brown dwarf companions have 
the same formation mechanism (see, e.g., Kouwenhoven et 
al. 2007a for a discussion). There are three methods of treat- 
ing these low-mass companions: 

- Accept all companions (PCP-I/SCP-I). All companions 
are accepted, irrespective of their mass. In this case, even 
stars with a planetary companion are considered to be in a 
binary system. 

- Reject low-mass companions (PCP-II/SCP-II). All com- 
panion stars with a mass smaller than a certain minimum 
mass Afiniii are rejected, and the primary star becomes a 
single star For this pairing function, stars with a low-mass 
(planetary) companion are not considered as binary systems. 
Note that the latter manipulation decreases the binary frac- 
tion, and increases the average mass ratio, particularly among 
low-mass binaries. 



-Redraw low-mass companions (PCP-III/SCP-III). When 
a companion with M2 < Mmin is drawn, it is rejected, and 
a new companion star is drawn. This procedure is repeated 
until M2 > Mmin- This pairing mechanism ensures that 
stars with low-mass (planetary) companions are not present 
in the resulting population, and that the resulting binary 
fraction is equal to B. This manipulation increases the aver- 
age mass ratio, particularly among low-mass binaries. 

4 Differences between pairing functions 

In this section we describe the main differences between the 
pairing functions described in §[3] Most of the properties of 
these pairing functions can be derived analytically. We refer 
to Kouwenhoven et al. (2008) for full derivations. 

For pairing function RP, the stellar mass distribution (of 
primaries, companions, and single stars) fM,a\\{M), is equal 
to the generating mass distribution fniM). For all other 
pairing functions, fM,a.\i{M) contains significantly more low- 
mass stars than fM{M). Only for SCP-1, the system mass 
distribution Jmi.m2.s{M) is equal to Jm{M). For pairing 
functions PCRP, PCP-1, PCP-II and PCP-III, the distribu- 
tion of primary/single stars masses (i.e., the target stars in a 
survey forbinarity) is equal to fmiM). 

The mass ratio distribution resulting from each pairing 
function is different. Moreover, for all pairing functions (ex- 
cept PCP-1) the mass ratio distribution varies with primary 
mass. This means that, for a given pairing function, the mea- 
sured mass ratio distribution and binary fraction depends on 
the range in spectral type of the selected targets in a binary 
survey. An example of the mass ratio distribution for the 
nine different pairing functions is shown in Fig.[Tl for three 
different targeted samples (all, high-mass, and low-mass tar- 
get stars). For the other pairing functions, the overall mass 
ratio distribution is biased towards low values. 

Each pairing function is characterised by a generating 
binary fraction B, the fraction of stars that is assigned a 
companion, prior to further manipulations associated with 
each pairing function. For most pairing functions the re- 
sulting overall binary fraction Bau is equal to B, except for 
PCP-11 and SCP-II, where Sail < B due to the rejection of 
low-mass (planetary) companions. The binary fraction is a 
function of primary mass for RP (except when B = 100%), 
PCP-11 and SCP-l/lI/111. Examples on the varying binary 
fraction for these pairing functions are shown in Fig.|2] 

5 Implications and interpretation 

The outcome of star formation theories can be described 
using pairing functions. Disk fragmentation, for example, 
may result in a pairing function similar to PCP-11, while the 
fragmentation of pre-stellar cores may result in a population 
similar to that of SCP-11. After (and during) the star for- 
mation process, the pairing function may be altered due to 
gravitational interactions, as well as stellar evolution, which 
can alter the properties of binary systems. 



© 2008 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-joumal.org 



Astron. Nachr. / AN (2008) 



3 



PCP-I PCP-II PCP-III RP PCRP SCP-I SCP-II SCP-III 



























L 












C2 













0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 

Mass ratio Mass ratio lyass ratio Mass ratio Mass ratio Mass ratio Mass ratio i^ass ratio 



Fig. 1 The mass ratio distributions resulting from the different pairing functions, for models with a Kroupa (2001) 
generating mass distribution, where applicable a generating mass ratio distribution fq{q) = 1, and a generating binary 
fraction of 100%. From top to bottom, the panels show the overall mass ratio distribution (i.e., for all binaries), the mass 
ratio distribution for binaries with 1.5 Mq < Mi < 20 Mq, and for binaries with 0.02 M© < AIi < 0.08 M©. This figure 
illustrates that each pairing function results in a different overall and specific mass ratio distribution. Note that, apart from 
the differences in the mass ratio distributions, there are also differences in the resulting mass distribution and overall and 
specific binary fractions (see Fig.|2]i. 



The easiest method to distinguish observationally be- 
tween the different pairing functions is to measure the mass 
ratio distribution as a function of primary spectral type. Us- 
ing this method, pairing functions RP and PCRP have been 
ruled out by observations (see, e.g., Kouwenhoven et al 2008, 
and references therein). A realistic pairing function should 
be able to account for the large number of massive, equal- 
mass binaries ("twin binaries") among massive stars (e.g., 
Pinsonneault & Stanek 2006; Lucy 2006; Soderhjelm 2007), 
as well as the small number of brown dwarf companions 
among stellar-type objects (e.g., Kouwenhoven et al. 2007a). 

Surveys for binarity have indicated that the binary frac- 
tion depends strongly on spectral type (see, e.g., Sterzik & 
Durisen 2004; Kohler et al. 2006; Lada 2006; Bouy et al. 
2006). For early-type stars (O/B/A) the binary fraction is 
close to 100%. The binary fraction decreases to 50 — 60% 
for F- and G-type stars, and decreases further to 30 — 40% 
for M-type stars. The binary fraction is lowest among brown 
dwarfs (10 — 20%). Assuming that this trend is not the result 
of selection effects, this trend is inconsistent with PCP-I and 
PCP-III (for which B is independent of Mi) nor SCP-I and 
SCP-III (for which B only mildly increases with Mi). Fur- 
thermore, observations of the mass ratio distribution have 
ruled out RP and PCRP (see above) in most stellar pop- 
ulations. Of the pairing functions described in these pro- 
ceedings, only PCP-II and SCP-II remain viable options. A 
deeper analysis, including a study of more complicated pair- 
ing functions, combined with further deep observations, is 
necessary for a full description of the pairing function in the 
different stellar populations. 

When deriving a pairing function from observations, a 
complication is introduced by selection effects. The origin 
of these selection effects is two-fold: sample biases (intro- 



duced by the choice of the observed sample), and instru- 
mental biases (e.g., due to limited spatial resolution or sen- 
sitivity); see Kouwenhoven (2006) for details. The instru- 
mental biases often prevent the detection chances of low- 
mass companion stars in a survey for binarity. Such selec- 
tion effects are usually well understood and can be taken 
into account. The selection of the targets in a binarity sur- 
vey introduces a sample bias. In the sections above we have 
seen that the mass ratio distribution and binary fraction de- 
pend strongly on the choice of the sample. In order to prop- 
erly constrain the pairing function of a stellar population, 
a survey of stars with a large range in spectral types is a 
prerequisite. To obtain an estimate of the primordial pairing 
function, one additionally needs to correct for the effects of 
dynamical and stellar evolution that have altered the binary 
population over time. 

6 Summary and discussion 

We have described several simple algorithms of pairing in- 
dividual stars into binary systems, and refer to these al- 
gorithms as pairing functions (Kouwenhoven et al. 2008). 
Each pairing function is characterized by a generating mass 
distribution fM{M) and a generating binary fraction B, and 
most pairing functions by a generating mass ratio distribu- 
tion fq{q). Several pairing functions can result in a binary 
fraction that differs significantly from the generating binary 
fraction B. The major differences are the dependence of the 
specific mass ratio distribution and the specific binary frac- 
tion on the selected primary mass range (see Figs[T]and|2]i. 

The pairing functions in these proceedings are described 
because of their simplicity or frequent use in literature. The 
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Fig. 2 The specific binary fraction Bmi{Mi) as a function of primary mass for pairing functions RP (left), PCP-II 
(middle) and SCP-II (right). Models with a generating binary fraction B of (from bottom to top) 5%, 30%, 60%, and 90% 
are shown. The minimum companion mass is set to 0.02 Mq. The horizontal lines represent the generating binary fraction. 
No selection effects have been applied. Left: results for RP with the Salpeter mass distribution. Middle and right: results 
for PCP-II and SCP-II with mass ratio distributions fq{q) oc q'*, with 7 = (solid curves), 7 = —0.3 (dash-dotted curves), 
and 7 — —0.6 (dashed curves). The figure shows that, even though we have applied the pairing function regardless of the 
primary spectral type, the binary fraction may vary strongly with primary spectral type. 



pairing function in real stellar populations could be differ- 
ent. The pairing function resulting from the star formation 
process may have an environmental dependence, and is al- 
tered over time by dynamical interactions between stars. 
Kroupa (1995a, 1995b, 1995c) and Thies & Kroupa (2007), 
for example, find that young binary populations are well de- 
scribed as a mixture of two different binary populations, 
both of which are randomly paired from restricted mass 
ranges. Note that such a heterogeneous population result- 
ing from restricted random pairing (RRP) over limited mass 
ranges, is significantly different from the population that re- 
sults from random pairing over the full mass range (RP). 

Each pairing function described in this paper results in a 
different mass ratio distribution and binary fraction. These 
properties also depend strongly on the targeted mass range 
in a survey for binarity. In these proceedings we have shown 
that it is important to carefully study selection effects in ob- 
servations, and to clearly state the pairing mechanism used 
in computer simulations, in order to draw conclusions on the 
overall binary population. The pairing functions described 
above are used in literature, but may not accurately describe 
realistic stellar populations. The next step towards under- 
standing the star formation process is to fully characterize 
the binary population in young stellar groupings. Variations 
in the primordial pairing function (e.g., binary fraction or 
mass ratio distribution) with environment necessarily imply 
variations in the IMF and star formation mechanism. 
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